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Abstract. When submitted to gentle mechanical taps a granular packing slowly compacts until it reaches 
a stationary state that depends on the tap characteristics. The properties of such stationary states are 
experimentally investigated. The influence of the initial state, taps properties and tapping protocol are 
studied. The compactivity of the packings is determinated. Our results strongly support the idea that the 
stationary states are genuine thermodynamic states. 
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PACS. 45.70.-n Granular systems, 
materials; granular materials 



45.70.Cc Static sandpiles; granular compaction - 81.05.Rm Porous 



1 Introduction 

In the absence of an external drive, granular materials 
rapidly come to rest. This is a consequence of their dis- 
sipative interactions and of the irrelevance of thermal en- 
ergy (in this case thermal energy is negligible compared 
to the energy needed to move a grain). Thus standard 
thermodynamics is not applicable to those systems. Ed- 

, wards and Oakeshott [T] proposed a non-standard thermo- 
dynamic description for static granular media. This theory 
postulates that a granular packing at rest can be described 
by suitable ensemble averages over its blocked "jammed" 
states. If one assumes that all the available blocked states 
have the same probability to occur and that the volume 

■ is analog to the energy of thermal systems, the configura- 
tional entropy of a granular packing is S = AlnJ7(y, N), 
where A is the equivalent of the Boltzmann constant and 
0(V, N) the number of mechanically stable configurations 
of N particles in the volume V. A temperature-like state 
variable, the compactivity X — dV/dS can then be in- 
troduced. This theory was partially investigated by recent 
experiments [2,3,4.5 , numerical simulations of models [6J 
I71I51IMTU] . These experiments have established that a gran- 
ular system subjected to a tapping dynamics slowly com- 
pacts (i.e. the packing fraction increases) and reaches a 
stationary state that depends on the tapping intensity. 
The existence of a stationary state is the first essential 
step to justify the validity of a thermodynamic-like de- 
scription for granular packings. 

Here we experimentally study the dynamics of granular 
media submitted to gentle mechanical taps. We mainly 
focus on the properties of the stationary states and dis- 
cuss the validity of such configurations as genuine ther- 
modynamic states. Using the fluctuations of the packing 
fraction at stationary state we will determine the com- 



pactivity introduced by Edwards and Oakeshott [T] and 
compare our results with those obtained for granular pack- 
ings submitted to fluid flow pulses [TT] . 
The outline of this paper is the following. We first describe 
the experimental setup. Then, in section [3j we rapidly re- 
call the previous results obtained on the relaxation of the 
packing fraction during compaction. Section 3] is devoted 
to the influence of the initial conditions on the stationary 
state and section [5] reports the dependence of our results 
on the frequency used to generate the taps. The cases of 
annealing and memory effects are studied in section [6] Fi- 
nally, section [7] is devoted to the conclusions of this work. 



2 Experimental setup 

The experiments (Fig. [It) are performed with d — 1 mm 
diameter glass spheres placed in a glass cylinder of di- 
ameter D = 10 cm. The cylinder containing the grains is 
tapped vertically at regular intervals (At = I s). Each tap 
is controlled by an entire cycle of a sine wave at a fixed 
frequency / = 30 Hz: V(t) = V M ax{1 - cos {2nft))/2 for 
< t < l/f and V(t) = elsewhere (Fig. [Tb). This ap- 
plied voltage is connected to an electromagnetic exciter 
(LDS V404) which induces a vertical displacement to a 
moving part supporting the container and the beads. The 
resulting motion of the whole system is monitored by an 
accelerometer at the bottom of the container. This mo- 
tion is more complicated than a simple sine wave: at first 
the system undergoes a positive acceleration followed by 
a negative acceleration with a minimum equal to —Jmax ■ 
After the applied voltage stops, the system relaxes to its 
normal position (Fig. [it). When ^ max is large enough, the 
bead packing takes off from the bottom of the container 
and achieves a flight until it crashes back to the bottom. 
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This crash is visible on the signal of the accelerometer: 
a negative acceleration — 7/ a z; corresponding to the fall of 
the grains is clearly visible. The control parameters are the 
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In a pioneering paper, Knight et al. [2] in Chicago 
first considered packing fraction relaxation law in gran- 
ular compaction. They showed that, starting from a loose 
packing of beads confined in a very narrow and tall tube 
(diameter 1.88 cm and 87 cm high, bead diameter 2 mm), 
a succession of vertical taps induces a progressive and very 
slow compaction of the system. They found that the re- 
laxation law can be well fitted by an inverse-logarithmic 
law, the so-called Chicago fit: 
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Fig. 1. a) Picture of the experimental setup, b) the signal 
controlling the tap (here f — 30 Hz) and c) the acceleration 
monitored by the accelerometer (here r » 3.3). 



frequency / used to generate the tap and the tapping in- 
tensity r — 'jmax/g, where g = 9.81 m.s -2 . The packing 
fraction is measured using a 7-ray absorption setup [3]. 
The measure is deduced from the transmission ratio of 
the horizontal collimated 7 beam through the packing: 
T = A/Ao, where A and Aq are, respectively, the activities 
counted on the detector with and without the presence of 
the granular medium. From Beer-Lambert law for absorp- 
tion, we can derive an estimation of the volume fraction 
in the probe zone: p « — (p,D)~ 1 fn(T), where p is the 
absorption coefficient of the beads. It was evaluated ex- 
perimentally to p, » 0.188 cm -1 for our 7 beam of energy 
662 keV ( 137 Cs source). The collimated 7 beam is nearly 
cylindrical with a diameter of 10 mm and intercepts per- 
pendicularly the vertical axis of the vessel. In order to 
reduce the relative fluctuations due to number of emitted 
photons we use an acquisition time of 300 seconds for each 
measure which leads to a precision of 0.001. The packing 
fraction of the sample is estimated from the ratio T av- 
eraged on approximately 7 cm height from the bottom of 
the cylinder. 



3 Packing fraction relaxation laws 

The first quantity of interest in compaction is the packing 
fraction (or density) p, defined as the ratio of the vol- 
ume of the grains to the total volume occupied by the 
packing. A few characteristic values of p for mono-sized 
sphere packings have to be reminded. The maximal pack- 
ing fraction reached in a random packing of spheres (the 
so-called random close packing fraction) is prop ~ 0.64. 
This value is significantly lower than the maximal pack- 
ing fraction obtained for face-centered-cubic (or hexago- 
nal compact) packing (p ma x ~ 0.74). Another limit is the 
so-called random loose packing corresponding to the less 
dense mechanically stable packing (prlp ~ 0.58). 
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For a given frequency /, the fitting parameters poo, p$, 
B and r essentially depend on r. The small number of 
grains in a tube diameter (ps 10) allows for a local mea- 
surement of the packing fraction with a capacitive method 
and prevents any convection in the packing. Nevertheless 
it induces strong boundary effects that may be responsible 
for crystallization (some packing fraction values obtained 
are well above the random close packing limit) visible in 
some of the experiments reported in 12 and in the values 
of poo obtained in [2] . 

More recently Philippe and Bideau [3] carried out new 
compaction experiments using the setup described in sec- 
tion [2] with about 100 grains in the tube diameter. This 
restricts the boundary effects but, contrary to the Chicago 
group's experiments, allows convection. The relaxation law 
obtained by these authors differs significantly from those 
obtained by Knight et al. [2], especially for the long-time 
behavior. Indeed whereas in previous experiments no clear 
evidence of convergence to a steady state has been estab- 
lished, such an evidence in definitely produced by our ex- 
periments and may correspond to a dynamical balance be- 
tween convection and compaction. The relaxation is better 
fitted by the Kohlrausch Williams Watts law (KWW law) 
- a stretched exponential - : 

p(t) = Poo~ {Poo - Po) exp [~(t/r f] 

where poo and po correspond respectively to the steady 
state and to the initial packing fraction value (figure [2]) . 
The adjustable parameters r and j3 are here respectively 
the relaxation time and a parameter related to the stretch- 
ing of the exponential. This characteristic time scale is 
found to be well described by an Arrhenius behavior 
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Such kind of relaxation law is also found for strong glasses 
(the dimensionless acceleration r plays the role of the 
temperature). Note that, Lumay and Vandewallc recently 
recover KWW-like laws for compaction of 2D granular 
packings 13.]. As can be seen in figure a stationary 
state is reached (if the number of taps is large enough). 
Figure [3] reports the evolution of the packing fraction of 
the stationary state as a function of the tapping intensity 
F for / = 30 Hz. We observe a decrease of the packing 
fraction with r. To reach an important value of packing 
fraction it is thus essential to use a small tapping intensity. 
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Fig. 2. Packing fraction versus the number of taps in our 
experiments. The curves are KWW fits (see text). The tap 0.6, 
frequency is / = 30 Hz. 
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Fig. 3. Packing fraction of the stationary state versus the 
tapping intensity for / = 30 Hz 



4 Influence on the initial state 

The question of the existence of a stationary state dur- 
ing granular compaction is essential to validate any ther- 
modynamic approach. As reported in the previous sec- 
tion, loose packings submitted to gentle mechanical taps 
continuously compact until they reach a stationary state. 
Does this state depend on the initial conditions ? To an- 
swer this question, we submit two granular packings of 
different packing fractions (one dense with a packing frac- 
tion close to that of the random close packing : packing 
1 and one loose with a packing fraction close to that of 
the random loose packing : packing 2) to taps of the same 
intensity. As reported in figure [4^ (see also [14] ) after a 
transient the two packings reach the same value of packing 
fraction. This behavior is observed for a large range of F 
(figure [4)3). This shows that, for a given tapping process, 
the value of the packing fraction at stationary state does 
not depend on the initial conditions. For the two initial 
packings considered, the number of taps needed to reach 



Fig. 4. a) Evolution of the packing fraction versus the num- 
ber of taps for an initially dense packing and an initially loose 
packing. We observe that the final packing fraction value is 
similar for the two packings (F — 2.1 and / = 30 Hz), b) Evo- 
lution of the final packing fraction versus the tapping intensity 
for an initially dense packing and an initially loose packing. 



the steady state is larger for decompaction than for com- 
paction. This is true for any value of r [14]. This can be 
intuitively understood : contrary to decompaction, com- 
paction is facilitated by the gravity. 



5 Influence of the tap duration 

As mentioned in previous section, the stationary state 
does not depend on the initial conditions. We can modify 
the frequency generating the tap and study its effect on 
the stationary state. We report in Fig the stationary 
values of the packing fraction p^ as well as its fluctua- 
tions Apao recorded after a sequence of taps of duration 
1/f (where / is 30 Hz, 60 Hz or 90 Hz) and intensity F. 



Even though p^ and its fluctuations at stationary state 
depend more or less on both r and 1/f (for the latter it 
is particularly visible at low values of r), when Apoo is 
plotted as a function of p^ the data collapse, within errors 
bars, onto a single master curve (figure EJd). This shows 
that poo and zipoo are linked. Note that the smallest values 
of Apao might be perturbated by the fluctuations of the 7- 
ray source. In other words, no matter how the state with 
packing fraction p^ is attained. This result is in agree- 
ment with those obtained very recently by Pica Ciamarra 
et al. [15] with numerical simulations of granular packing 
submitted to fluid flow pulses. These results strongly sup- 
port the idea that such stationary state are indeed genuine 
"thermodynamic states". Moreover, in [T5], the authors 
demonstrate that at stationarity a granular packing can 
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order polynomial fit (figure [6| . In order to compare our 
results with those reported in [11) we use the same value 
for the packing fraction of the random loose packing : 
Prlp = 0.573. We observe a decrease of X with p^. al- 
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Fig. 5. a) The main panel shows the packing fraction obtained 
at stationarity as a function of the tapping acceleration for 
different values of the frequency of the mechanical tapping. The 
inset shows the dependence of the packing fraction fluctuations 
on r and /. b) Plot of the packing fraction fluctuations versus 
the packing fraction at stationarity. All the data collapse on a 
single master curve, showing that Ap x and p x are linked and 
that this link does not depend on the values of / and F. 



be described by only one parameter, the packing fraction, 
all the other observables being characterized by this pa- 
rameter. The strong resemblance of our results with the 
ones reported in [15] suggests that this is also true for ex- 
periments of granular packings submitted to mechanical 
tapping. 

It should be pointed out that this is true only at the sta- 
tionary state. Indeed, during the transient, memory effects 
can be observed [TB]. In this case, the packing fraction is 
in general no more enough to fully characterize the pack- 
ing and other observables are needed [17]. Knowing the 
dependency of Ap^ on p^ enables us to determine Ed- 
wards' compactivity, which the equivalent of the tempera- 
ture in thermal systems. Using the fluctuation dissipation 
theorem 1181 we can derive 



^X(p) = 
m 



' PRLP 



where we have used X(prlp) = oo. In this expression 
m is the grain mass and p g the grain density. We an- 
alytically solve this equation using for Ap 00 (p 00 ) a 2 nd 




Fig. 6. Compactivity as a function of the average packing frac- 
tion at stationary state. (Inset) Evolution of the compactivity 
versus the tapping intensity for / = 30 Hz. 



though the evolution of Apoo is different, the behavior of 
X(p oc ) is similar to that found for granular packings sub- 
mitted to fluid flow tapping [11] . We also report in the in- 
set of figure [6] the evolution of the the compactivity versus 
the tapping intensity r for / = 30 Hz. We observed that 
these two quantities are linked: the temperature- like com- 
pactivity increases with the tapping intensity. Note that 
the change of behavior observed around r — 2 may be ex- 
plained by the two different convective regimes observed 
in our system [3]. In the analogy between glassy systems 
and granular packings undergoing compaction, r is sup- 
posed to play the role of temperature at the stationary 
state. So the above-mentioned result reinforce this anal- 
ogy. It should be interesting to compare the compactivity 
with other "temperatures" defined for granular media as 
the granular temperature (defined as the velocity fluctu- 
ations) or the effective temperature defined through the 
out-of-equilibrium fluctuation-dissipation theorem [19 . 



6 Annealing and memory effects 

Further insights into the understanding of the stationary 
state characteristics can be gained by allowing the tap 
intensity to vary in time. Annealing experiments phenom- 
ena in granular media undergoing compaction have been 
reported previously in [12] . The authors proceeded as fol- 
lows: starting from a loose packing of grains, the mate- 
rial is tapped at a given intensity r for a given time t 
(10 5 taps), r is then modified and the compaction pro- 
cess continued for t taps (see figure [7^,) . The increase of 
r corresponds to an increase of the average packing frac- 
tion except for values larger than three for which a slow 
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1 we can study the response of that packing to taps char- 
acterized by F 2 and f 2 and compare with the behavior 
of the packing at stationary state in experiment 2. Since 
the history of the two packings are different one may ex- 
pect a transient in the behavior of the packing 1, sign 
of memory effects. Using this protocol no packing frac- 
tion modification can be seen. So no memory effects exist 
for packings at stationary state. The range of parameters 
explored remains small but this preliminary result con- 
firms that the packing fraction is an observable that fully 
characterize the stationary state. Once more, this not true 
during transient when memory effects |16j reflect the need 
for, at least, an extra-observable to fully characterize the 
packing. 







12 



Fig. 7. a) Sketch of the evolution of the packing fraction during 
annealing experiments reported in [12] . b) Results of annealing 
experiments in our experiments with / = 30 Hz. No aging is 
observed. 



decrease can be observed. If r is then reduced, the pack- 
ing fraction, rather than following reversely the previous 
curve, still increases. This new curve, contrary to the pre- 
vious is reversible. Using our setup and the same protocol 
(with / = 30 Hz), we only recover the reversible branch 
(see figure [7b). This can be explained by the fact that the 
number of taps used is large enough to allow our system to 
reach stationarity. Thus, aging and irreversible-reversible 
behaviors are observed only when the stationary state is 
not reached. Note that these experiments also confirm that 
the stationary state only depends on the tapping protocol 
and not on the initial conditions. 

We also carried out memory effects experiments. Let 
us recall that Josserand et al. [16] carried out experiments 
on the response function of a granular media undergoing 
compaction to a sudden perturbation of the tapping inten- 
sity. These authors drive a granular packing to the same 
packing fraction po with three different tapping intensi- 
ties. Then those three packings are tapped at the same 
intensity Jg and their behavior depends, for short time, 
on the previous value of the tapping intensity. These data 
show a short-time memory effect; the future evolution of 
the packing fraction depends not only on its initial value 
but also on its history. Using our experimental setup we 
recover those memory effects as reported in [20l[T7]. Do 
such memory effects also exist at the stationary state ? 
We can carry out two compaction experiments until the 
stationary state is reached. For the two experiments we 
use two different tapping characteristics (i.e. frequency /i 
and tapping amplitude i~i for experiment 1 and frequency 
f 2 and tapping amplitude F 2 for experiment) in such away 
that the packing fractions of the two stationary states are 
equal. Once the stationary state is reached in experiment 



7 Conclusions 

To Summarize, in this paper we show that granular pack- 
ings submitted to gentle mechanical taps can reach a sta- 
tionary configuration. This state does not depend on the 
initial conditions. A dependence on the duration of the 
taps is found but our result show a one-to-one correspon- 
dence between the final packing fraction and the packing 
fraction fluctuations. The compactivity, a temperature- 
like state variable, is then determinated as a function of 
packing fraction of the stationary state and as a function 
of the tapping intensity at a given /. We also show that 
neither memory effects nor annealing are present at sta- 
tionary state. Our results strongly support the relevance 
of a fundamental theory of dense granular media. 
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